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^ ; Abstract 
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■ Various scenarios of the initial inflation of the universe are distinguished 

■ by the choice of a scalar field potential U{(j)) which simulates a temporar- 

u '. 



ily non-vanishing cosmological term. Our new method, which involves a 
reparametrization in terms of the Hubble expansion parameter H, provides a 



^ ■ classification of allowed inflationary potentials and of the stability of the crit- 



ical points. It is broad enough to embody all known exact solutions involving 
one scalar field as special cases. Inflation corresponds to the evolution of criti- 
cal points of some catastrophe manifold. The coalescence of its nondegenerate 
critical points with the creation of a degenerate critical point corresponds the 
reheating phase of the universe. This is illustrated by several examples. 
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Introduction. In the 80 's, Guth |T[ and Linde have modelled an inflationary phase 
of the universe (cf. 0). Scalar fields (Higgs, axion) are expected to generate, shortly after 
the big bang, an exponential increase of the universe. The so-called graceful exit to the 
Friedmann cosmos was partly solved in the new inflationary universe 0. In this model, 
the scalar field is ruled by a slightly different self-interaction potential which possesses a 
slow-roll part (a plateau) of the potential (acting as a vacuum energy) which dominates the 
universe at the beginning. Power-law models were constructed which possess no exponential 
but an a{t) ~ increase of the expansion factor of the universe The intermediate in- 

flation is merely a combination of exponential and power-law increase Further solutions 
were found in 0. Recently, by using the Hubble expansion parameter if as a new "time" 
coordinate, we [|1^ were able to derive the general Robertson-Walker metric for a spatially 
flat cosmos. Our formal solution for arbitrary U (0) comprises all previous exact solutions. 

For a rather general class of inflationary models the Lagrangian density reads 
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/: = —yJ\g\\^R + K[g^''{d,<P){d,<P) - 2f/(0)J j , (1) 

where is the scalar field and U{(f)) the self-interaction potential. We use natural units with 
c = h = 1. A constant potential Uq = A/k, would simulate the cosmological constant A. 



We do not separately discuss non-minimally coupled Jordan-Brans-Dicke type models ||Tl 
since they can be reduced to (|l|) via the Wagoner-Bekenstein-Starobinsky transformation 
r^^lSl . Let us concentrate on the flat Friedmann-Robert son- Walker cosmos 



ds^ = dt^ - a\t) [dr^ + (dO^ + sin^ Odip^)] , (2) 

where a{t) is the time-dependent expansion factor with the dimension length. Flat space is 
anyhow favored in the inflationary scenario. The scalar field depends only on the time t, 
i.e. (p = (pit). 

The Friedmann evolution equations. Let us assume that a{t) ^ 0, such that we 



can express our equations completely ||T0[ in terms of the Hubble expansion rate if := 
a{t)/a{t). Only the diagonal components of the Einstein equation are non-vanishing. The 
(0, 0) component involving the density p, reads 



3H'^ = Kp = + U] . 



(3) 



It describes the conservation of the energy. The (1,1), (2,2), and (3,3) components are 
given by 

'1 V 



2H + 3H^ = -Kp = -K{-(j)'^ -U 



(4) 



where p is the pressure generated by the scalar field. The resulting Klein-Gordon equation 



IS 



-3H<j) - U'{(j)) 



(5) 



which is, after multiplication by 4>, 

By linear combination of (|^) and (^ we find the autonomous nonlinear system 



(6) 



H = KU{(j)) - 3H^ =: V{H, (f)) , 
= ±J'^^j3H^ - KU{(t)) = ±^ 



K 



'-V{H, 



(7) 
(8) 



The function V{H,(f)) will turn to be the "height function" in Morse theory [T^. Observe 
that (H) is, in view of (^ and (^, a first integral of (§). Moreover, < in order to avoid 
scalar ghosts. For the metric (0), the Lagrangian density (P reduces 



C 



(9) 



Since the shift function is normalized to one for the metric (0), the canonical momenta are 
given by p = dC/da = —GHo^/k and vr = dC/dcj) = a? (p. [Incidentally, this suggest to 
take the volume as a generalized coordinate and use P = dC/d{a^) = —2H/k, as new 
momentum.] The Hamiltonian or "energy function" is given by 



E 



(10) 
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and vanishes for all solutions. 

Catastrophe of Whitney manifold and the critical points of the evolution. 

In the phase space [0, the equilibrium states of the system (|^ and (H) are given by the 



constraint {H, 0} = 0. The critical or equilibrium points, respectively, of this system are 
determined by V{Hc, 0c) = 0. This constraint is globally fulfilled by K,U{(f)) = 3H^, where 



the Hubble expansion rate is constant, i.e. Ha ='■ y A/3 . For = and A 7^ 0, we obtain 



the de Sitter inflation with a{t) = aoexp(yA/3t). 

The Jacobi matrix J of the system (0) and (||) is given by 

/ -6H kU' \ 

J={ h (11) 

V±6i/(-2/tV)-i/2 ^/t[/'(-2/tV)-i/2y 

where U' = dU/dcj). Since det J = 0, the system is degenerate. For the analysis of stability, 
it suffices therefore to consider only (|^) and later to reconstruct 0. Since H and are 
independent variables, we can introduce the non-Morse potential W{H, (p), which is defined 
via V := —dW/dH and analyse the system with the aid of catastrophe theory. From (|^ we 
obtain 

W{H,^)=H''-kU{(P)H + C{(P) , (12) 

where C is an arbitrary function of 0. The function W is already in canonical form in 
if-space, and belongs to a Whitney surface or to the Arnold singularity class A2 (cf. [0]). 
That means that our catastrophe manifold here is the Whitney surface. This manifold has 
only one control parameter, which is here given explicitly as the potential U. Thus, an 
evolution of critical points are determined via the values of the potential U. Let us analyse 
the types of critical points at different fixed values of the control parameter U. 

If Uc '■= U{(j)c) < 0, the equation has no stable critical point due to the shape of the 



Whitney surface. However if Uc > 0, there are two critical points: stable at He = JUd^ 



and unstable at He = —yUc/<i. For He > 0, Walliser ||T9| comes to the same conclusion. 

The value ?7c = is the bifurcation point. Provided this is also an extrema of V, we 
necessarily have dV/dH |c= —QH^ = 0, dV/dcf) \c= kU'^ = 0, and 0c = 0. Thus, also the 



critical points of the Klein-Gordon equation are involved. Hence, the Hubble parameter has 
to vanish and (pc is a double zero of the potential U. The Hessian of (|^) takes the form 

f-6 \ 

Hess{V) = \ . (13) 

V kU"J 

The sub-determinant of the Hessian is Aq = d^V/dH"^ = —6 < and Ai = det Hess{V) = 
(d^V/dH^) {d'^V/d(f)'^) - {d'^V/{dHd(f))f = -6kU". For a maximum of the potential U, 
i.e. U' = and U" < 0, the function V possesses a maximum; for a minimum of the 
potential U we find a saddle point for V. 

Alternatively, we can investigate the non-Morse potential (|l^) alone, e.g. for the chaotic 
inflationary model with U{(f)) = 0^/ n 

W{H,(P)=H''~4>^H + C{cP) , (14) 

which has minimum and maximum at He = ±(f)c/V^ and a saddle point for 0c = 0. This 



corresponds to the end of inflation and a reheating of the universe pO . 

In the case of power-law inflation, we have U (0) = (exp 4>)/k and therefore obtain 

W{H,(j)) = -e'f'H + C{(j)) , (15) 

which has the critical points = ±e^''^'^/\/3 and a saddle point for 0c = —oo. The latter 
can be reached because of = ln(l/t) for t —>■ oo. In both models, the saddle point appears 
at that time of the reheating of the universe. 

The last example here should be taken from the new inflationary model with t/(0) = 
[0^(0^ — A) + D]/k,, where A, D are two constants. The function 

W{H, (p) = H^- [0'(0' -A)+ D]H + C(0) (16) 



possesses minima and maxima at He = ±yU{(f)c)/S and two saddle points at 0c = ±y A/2, 
the zeroes of the potential, if D = A'^/A. That means that a shift of the potential is necessary, 
such that the zeroes of the potential are also its minima. 

The reheating phase of the universe. Since ?7(0c) corresponds to the latent heat of 
the universe in this phase, we are now in the position to state the following: 



Theorem: The critical points of the non-Morse potential W{H, 0) determine the evo- 
lution in the inflationary phase. Along the minima and maxima = ±-y/f/(0j/3, the 
inflaton moves from the slow-roll to the hot regime. The saddle points of W , i.e. more 
precisely, the minima of V, determine the onset of reheating. 

The critical points are related to some turning points of the evolutions. Since our po- 
tential U or the field 0, respectively, depend on time, we obtain the following scenario: The 
system starts from some initial conditions. In the first stage of the infiation, the stable criti- 
cal points define the regime of the evolution, that is the system evolves towards such points. 
In the second stage, because our control parameters (here the potential U or the field 0, 
respectively) also depend on time, the evolution moves along the critical points. Consider, 
for example, the chaotic infiationary model in Eq. (p!4D: The system being initially located in 
a minimum, evolves to the saddle point; physically, this means the end of the infiation and 
a reheating of the universe. In this way, we have derived a universality class for infiationary 
models. This picture is very general, as we are going to show in the next section by means 
of catastrophe theory. 

Bifurcation and other catastrophes. For the analysis of critical points of multidi- 
mensional functions (see |]rB|) the Hessian of (|T2|) is usually employed. It has the form: 

6H -kU' \ 

Hess{W) =1 . (17) 

-kU' -HkU" + C" 



Thus the classification of all critical points may be given by the Whitney theorem [18 



However, the explicit form of (|l^) allows to solve this problem without an analysis of the 
Hessian. The given form of W suggests that we have here umbilic catastrophes. There 
are only a few of them: elliptic or hyperbolic Df, second elliptic or second hyperbolic 
{Dq or Dq, respectively), and symbolic umbilic Eq. With the knowledge of the form of 



these catastrophes, which are given for example in Ref. [|T^] one can completely predict the 
evolution of the system. For the example of the chaotic infiation, Eq. ( [T^ ) describes an 
elliptic umbilic catastrophe which has, in this special form, only one critical point. 

For K,U{(f)) 7^ 3H^, we find {H, 0} ^ 0, which implies that the solutions = 0(t) and 



H = Hit) are invertihle, i.e. t = t(H) = f -. In contrast to the construction of 



Lidsey p2[ in which the scalar field itself is employed as a new time variable, our approach 
is also valid for the end of infiation. Then we can write the potential in (1^) and (M) in the 



reparametrized form [|T0 



u{<f>) = um)) = uimm = u{h) . (is) 

The reduced problem is given by the one-dimensional equation 

H = kU{H) - 3H^ = g{H) , (19) 



where g{H) is the "graceful exit function" of Ref. |jTO . 

All solutions of this equation may also be classified with the aid of the catastrophe theory 
||18|| . Eq. (|l^) has critical points which are determined via the function W{H), defined, 
analogously to W, by g := —dW/dH. The simplest case arises when W{H) is exactly the 
Morse potential, i.e. W{H) = XH'^/2, i.e. its critical points are nondegenerate. In this 
simple case the solution depends on the sign of A. If A > 0, the only possible stable state, 
to which the system evolves, is H = 0. 

For A < 0, we obtain H C exp \t and the system has no stable critical points. In that 
case a ~ Ci exp(cexp(At)). Then, the universe would be expanding too fast. 

The regimes of evolutions, or critical points, depend on the shape of the potential W{H). 
With the aid of the theory of singularities (in that particular case with the aid of elementary 
catastrophe theory |]l8l) we may classify the infiation regimes. If the potential W{H) is a 
smooth function it belongs to the one of the Arnold classes A^-, where n > 2. 

The canonical form of the n-th class is W{H) = XH"'^^. The co-dimension of that 
potential is equal to n — 1, that is, the number of control parameters is equal to {n — 1). 

The critical points, which are structurally stable correspond to the minima and maxima 
of the polynomial 

W{H),,form = A(if"+^ + Xn-lH''~' + A.-S^^""' + •.• + Aq) . (20) 
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The system described by ( |T9[) make an evolution to the minima of W{H) deform- In each of 
the minima, we have H = and H = const. The maxima of the potential (|20|) correspond 
to unstable points for which the instability transforms the system to its minima. With the 
change of the control parameters, some minima coalesce with maxima or vice versa. As a 
result, the number of minima changes. Each minimum is associated with a distinct constant 
value of the Hubble expansion rate. 

That is, each of the critical points corresponds to a different de Sitter type inflation. 
However, as in the considered case of the Whitney catastrophe the coalescence of these 
nondegenerate critical points into degenerate one mean the reheating phase of the universe. 
That is, locally, on the considered complicated manifold, we have a Whitney sub-manifold. 
The inflationary evolution and the picture of the creation of the reheating phase is the same 
as it is described above. 

However, there are a few exceptions related to the instabilities. One of these instabilities 
for the non-Morse potential we have already discussed above. Asymptotically, for a large 
time scale, the instabilities are defined by the leading term of the polynomial Wdeform{H). 

Asymptotically, if A < and n > 0, we may write for the n-th Arnold class 

g{H) = -{n + l)XH\ (21) 

such that ([T9|) yields the solution 

i/= [C+(n2-l)At]i/(^-") (22) 

where C is some positive constant. These instabilities have no physical meaning. 

Among the explicit models which have been analyzed in , we consider as an instructive 
example g{H) = —2H'^/A^ + 2A^\^ (in the notation of 0). This corresponds again to the 
Whitney catastrophe or to the second Arnold class A2. Therefore, we may expect the 
described above picture. Let us show it explicitely. This ansatz leads to 

0(t) = Aln[tanh(At)] , (23) 
H{t) = A^\cot\i{2\t) , (24) 
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a{t) = ao[sinh(2At)]^'/2 ^ (35) 



f/(0) = A^X^ 



(3A2-2)coshM^ +2 



(26) 



as solution. Thus we have recovered one of the recent solutions of Barrow ||9||, for further 



details see 10 



Now we have our two possible descriptions of the investigations of the critical points. 
The potential U (0) has to be shifted to 

U{(P) = A^X\3A^ -2)cosh^ , (27) 

so that the only critical point is at the origin = 0, which is a minimum of U and hence a 
saddle point of W. Alternatively, we can use g{H), which also has to be shifted into 

g{H) = -2H'^/A^ . (28) 

The non-Morse potential W = 2H^/3A^ possesses, too, a saddle point at = 0. In the 
two equivalent descriptions we were able to find that the saddle point corresponds to the 
reheating phase. 

Furthermore, it is now possible to relate the "slow-roll" condition, for the velocity of 
the inflationary phase, to the critical points resulting from catastrophe theory. For inflation 
(with a > 0) the two "slow-roll parameters" are given [|1^], in first order approximation, by 
e = —g/H^ and rj = —dg/d{H'^) , where g is the "graceful exit function" given in (|I^. In 



this reduced dynamics, they are effectively determined by the first and second derivatives of 
the reduced non-Morse function W{H), i.e., more precisely by e = —{l/H'^){dW/dH) and 
rj = {l/2H)d'^W/{dH)'^. They will also determine the density fluctuations |2l|. 
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